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THE  DIOCOTRON  INSTABILITY  IN  ANNULAR  RELATIVISTIC  ELECTRON  BEAMS 


Michael  E.  Jones  and  Michael  A.  Mostrom 

Intense  Particle  Beam  Theory  Group 
Los  Alamos  Scientific  Laboratory 
Los  Alamos,  New  Mexico  87545 

ABSTRACT 

The  instability  is  shown  to  have  a  finite  k  bandwidth  with  the 
fastest  growing  mode  in  the  lab  frame  shifted  to  nonzero  k.  The  thin 
beam  dispersion  relation,  the  appearance  of  several  new  parameters ,  and 
the  instability  mechanism  are  discussed. 


The  diocotron  instability  in  nonrelativistic  annular  electron 

beams  has  long  been  utilized  for  the  generation  of  microwaves  in  the 

1  2 

crossed  field  magnetron.  *  Allowed  to  evolve  to  a  nonlinear  state, 

3  4 

the  instability  can  disrupt  the  beam  by  causing  f ilamentation .  ’  The 

recent  interest  in  areas  employing  intense  annular  relativistic  elec- 

5  b 

tron  beams,  such  as  free  electron  lasers,  multi-gapped  accelerators, 

and  the  Anomalous  Intense  Driver  fusion  concept,^  makes  it  particularly 

important  to  understand  the  diocotron  instability  in  the  relativistic 

regime.  The  present  analysis  of  the  instability  provides  a  means  of 

determining  constraints  necessary  for  propagation  of  these  beams  down 

long  conducting  drift  tubes  and  for  evaluating  the  feasibility  of  using 

intense  relativistic  beams  for  high-power  microwave  generation. 

Previous  relativistic  analyses  of  the  diocotron  instability  have 

been  limited  to  the  special  cases  where  either  the  wave  frequency  u>  or 
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the  axial  wavenumber  k  have  been  set  equal  to  zero.  ’  A  complete 
relativistic  theory  must  of  necessity  include  U)  and  k  both  nonzero.  In 
the  analysis  presented  here  the  dispersion  relation  is  calculated  in 
the  rest  frame  of  the  beam.  The  dispersion  relation  in  the  lab  frame 
is  then  obtained  by  applying  a  Lorentz  transformation: 

=  vCu^+^v)  ,  k£  =  vC^+u^v/c2)  (1) 

where  the  £  and  b  subscripts  denote  quantities  evaluated  in  the  lab  and 
beam  reference  frames  respectively.  The  velocity  of  the  beam  is  v  and 
y=(l-v  /c  )  ^  where  c  is  the  speed  of  light.  From  Eqs .  (1)  it  is  clear 
that  in  order  to  obtain  the  lab  frame  quantities  it  is  necessary  to 
know  u)  as  a  function  of  k  in  the  beam  frame.  Therefore,  much  of  the 
present  analysis  concerns  the  extension  of  the  beam  frame  problem  to 
include  finite  k. 
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The  electron  beam  is  assumed  to  be  initially  in  an  azimuthally 
symmetric  equilibrium  inside  a  conducting  drift  tube  of  radius  d  as 
illustrated  in  Fig.  1(a).  A  large  externally  applied  axial  magnetic 
field  is  required  for  equilibrium.  In  the  beam  frame  the  only  equili¬ 
brium  motion  of  the  electrons  is  an  ExB  rotational  drift  induced  by  the 
applied  magnetic  field  and  the  radial  electric  self-field  of  the 
unneutralized  beam.  If  the  azimuthal  motion  is  nonrelativistic ,  then 
the  beam  frame  problem  reduces  to  the  case  considered  by  Levy.^ 

The  electrostatic  approximation  is  used  in  the  beam  frame.  This 
does  not  mean  that  the  perturbed  magnetic  field  is  zero  in  the  lab 
frame,  however,  since  the  Lorentz  transformation  of  the  electric  field 
of  the  beam  frame  will  give  rise  to  a  magnetic  component  in  the  lab 
frame.  Using  Maxwell 's  equations  and  the  linearized  cold  fluid  equa¬ 
tions  for  the  electrons,  solutions  of  the  form  E(r)exp{i(kz+£0-u>t) } 
yield  the  following  equation 

V  •  (e  •  E)  =  0  (2) 

where  V  denotes  the  operator  e^r  ^(d/dr)r  +  eg  i£/r  +  ik,  and  £  is 

the  dielectric  tensor  whose  only  nonzero  components  for  u)  <<  and 

2  2 

ui  «  U)  are  £  =  1,  £  Q  =  -£n  =  -iu)  /ww  ,  £Qn  -  1+tu  u)A/u)(£u)  -w)  x 

0  c  rr  9  r0  0r  p'  c>  00  p  0  0 

[(r/u>  )d(£n  U)2)/dr-2k2r2/£(£u)  -w)  ]  ,  and  £  =  l+ui^/ui(£u)n-u))'.  The 

c  p  u  zz  p  U 

plasma  frequency  of  the  beam  is  denoted  by  and  is  the  cyclotron 
frequency  associated  with  the  applied  magnetic  field.  The  equilibrium 
rotational  frequency  U)g  of  the  electrons  is  a  function  of  the  radial 
position  r.  Introducing  the  electrostatic  potential  4>  such  that 
E  =  ~V4>,  Eq.  (2)  becomes 
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r  *d(rd<|>/dr)/dr  -  (£2/r2  +  - 

cl>  is,  u)2  d  ?  k2tu2  ) 

j;  /  ^(*n  V  "  Uvi)|  (3) 

which  reduces  to  Levyfs  Eq.  (10)  for  k-0.1  A  similar  equation  has 
been  derived  in  a  study  of  an  instability  driven  by  shear  in  the  axial 
velocity.^ 

The  diocotron  dispersion  relation  is  obtained  by  solving  Eq.  (3) 
in  the  three  regions  shown  in  Fig.  1(a).  The  beam  density  is  nonzero 
only  in  region  II  (a  <  r  <  b).  The  potential  4>  is  required  to  be 
continuous  and  to  vanish  at  the  outer  conductor  and  at  the  origin. 
Jump  conditions  on  d4>/dr  are  obtained  by  integrating  the  differential 
equation  across  the  vacuum-beam  interfaces. 

Before  solving  for  the  complete  dispersion  relation,  it  is  in¬ 
structive  to  consider  the  physical  mechanism  of  the  instability.  For 
the  special  case  of  a  step-function  radial  density  profile  and  k=0 ,  the 
dispersion  relation  reduces  to  a  quadratic  equation  in  u)  corresponding 
to  two  modes.1  These  waves  cause  density  perturbations  only  at  the 
edges  of  the  beam  and  hence  are  surface  modes.  By  considering  the  even 
more  special  case  in  which  the  outer  surface  of  the  beam  is  in  contact 

with  the  conductor  (d=b),  one  finds  that  the  dispersion  relation  yields 

2Z 

one  mode  which  corresponds  to  a  wave  with  frequency  Q=l-(a/b)  propa¬ 
gating  on  the  inside  surface  of  the  beam.  Here  Q  is  the  frequency  10 
normalized  to  the  diocotron  frequency,  u)^=u)^/ 20)^ .  Similarly,  placing  a 
conductor  in  contact  with  the  inner  surface  of  the  beam  shorts  out 
perturbations  on  the  inner  surface  and  the  beam  again  has  only  one 
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2  2 

normal  mode  corresponding  to  a  wave  with  frequency  fi=£(l-a  /b  )- 

00  OQ  OO 

[l-(a/bj  ][l-(b/d)  ]/[l~(a/d)z  J  which  propagates  on  the  outside  sur¬ 

face  of  the  beam.  The  wave  energy  E  of  these  modes  is 

d 

0  . 

where  is  the  real  part  of  the  dielectric  tensor. From  this  expres¬ 
sion,  it  is  easy  to  show  that  the  inside  surface  wave  has  positive 
energy  and  the  outside  surface  wave  has  negative  energy.  The  disper¬ 
sion  diagram  for  these  modes  along  with  the  k=0  diocotron  dispersion 
diagram  is  shown  in  Fig.  1(b).  Thus,  it  becomes  clear  that  the  dioco¬ 
tron  instability  results  from  the  coupling  of  a  positive  energy  wave 
with  a  negative  energy  wave  in  a  manner  analogous  to  the  mechanism  of 
the  two-stream  instability.** 

It  has  been  shown  that  if  a  sufficiently  large  amount  of  charge* 
12 

or  current  (which  has  an  equivalent  effect  as  charge  in  the  beam 

frame)  is  placed  inside  the  annular  beam  the  diocotron  modes  can  be 

stabilized.  The  addition  of  a  charge  Q  increases  the  difference 

between  the  frequencies  of  the  outside  and  inside  surface  modes  by  an 
2  -2  -2 

amount  AQ=(4£Qe/mu)^)  (a  -b  ).  Thus,  if  these  modes  are  near  resonance 
for  Q=0,  then  the  addition  of  a  large  amount  of  either  positive  or 
negative  charge  results  in  a  frequency  shift  which  precludes  the  pos¬ 
sibility  of  the  coupling  of  the  positive  energy  wave  with  the  negative 
energy  wave. 

The  right-hand  side  of  Eq.  (3)  contains  resonance  terms  which  make 
ill e  differential  equation  singular  when  £u)  (r)=u).  By  considering  only 


/ 


dr  (E' 


•9(u£  )/9u> 


E} 


(4) 
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the  special  case  of  a  step  function  density  profile  and  k=0,  one  avoids 
the  question  of  these  resonances  since  the  coefficients  of  these  terms 
vanish.  However,  the  solution  for  this  case  reveals  that  the  diocotron 
mode  frequencies  are  precisely  in  the  right  range  of  values  to  be  reso¬ 
nant.  It  has  been  shown  that  for  a  smoothly  varying  density  profile 
the  first  term  of  the  right-hand  side  of  Eq.  (3)  can  lead  to  growth  or 
damping  depending  on  the  sign  of  the  term.^  The  physical  mechanism 
for  this  process  is  a  wave-particle  interaction  analogous  to  Landau 
damping.  The  second  term  on  the  right-hand  side  cannot  be  avoided  by 
choosing  a  specific  density  profile.  Because  this  term  is  always  nega¬ 
tive,  one  would  expect  from  the  analysis  of  the  contribution  of  the 
first  term  that  the  second  term  always  provides  a  damping  effect  for 
finite  k.  For  large  k  this  wave-particle  interaction  should  dominate 
the  wave-wave  interaction  and  the  instability  will  be  stabilized.  Thus, 
the  diocotron  instability  has  a  finite  bandwidth  in  k. 

An  upper  bound  on  the  growth  rate  and  bandwidth  may  be  obtained  by 
the  means  of  quadratic  forms. ^  Defining  t|>  by  4>=(£i0q-u))^  i|i,  multiply- 
ing  Eq.  (3)  by  \\i  and  integrating  from  r  =  0  to  r  =  d  one  obtains 

UK2  <  {(£u)Q/2)2-k2U)2}/{k2+(jA/d)2}  (5) 

where  ut  is  the  imaginary  part  of  U),  is  the  first  zero  of  the 

ordinary  Bessel  function  of  the  first  kind  and  order  £,  and  is  the 
derivative  of  U)^  evaluated  at  its  largest  value.  This  expression 
demonstrates  that  the  diocotron  modes  become  stabilized  if  the  shear  in 
the  equilibrium  rotation  frequency,  which  provides  the  free  energy  for 
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the  instability,  vanishes.  Equation  (5)  also  shows  that  the  instability 

must  stabilize  for  large  k  (the  expression  is  valid  only  for  UK  ^  0). 

For  the  step-function  beam  density  profile,  Eq.  (5)  restricts  the  range 

of  k  over  which  instability  can  occur  to  values  such  that  Ak  <  2H/v\ 

where  K=ka  and  n=2u)  /u)  . 

c'  p 

In  what  follows,  Eq.  (3)  is  analyzed  for  the  step-function  beam 
density  profile.  The  generalization  to  other  density  profiles  is 
straightforward.  It  should  be  pointed  out,  however,  that  other  pro¬ 
files  may  be  more  or  less  unstable  owing  to  the  Landau  damping  (or 
growth)  contribution  of  the  first  term  on  the  right-hand  side.  In 
regions  I  and  III  the  right-hand  side  of  Eq.  (3)  vanishes  and  the 
eigenfunctions  in  these  regions  are  expressed  in  terms  of  the  modified 
Bessel  functions  I^(x)  anc^  K^(x)  where  x^kr.  For  the  step-function 
density  profile,  w0=U)d(l-a2/r2) .  Thus,  in  region  II,  Eq.  (3)  has  a 
second-order  pole  at  r/a=p^=£/ (£-Q) .  Using  a  change  of  variables, 
u=£n(r/a),  Eq.  (3)  in  region  II  can  be  cast  in  the  form  of  a  Schroe- 
dinger  equation:  d  4>/du  -V(u)<J>=0.  Because  most  of  the  new  physics  for 
finite  k  results  from  the  resonance  term  or  pole,  V(u)  can  be  approxi¬ 
mated  by  expanding  about  this  pole.  This  expansion  should  be  particu¬ 
larly  good  for  a  thin  beam,  (b-a)/a<<l,  since  all  points  in  region  II 
are  then  close  to  the  pole.  Defining  a  new  variable  z=6  Jin(r/ap^), 
then  this  equation,  neglecting  terms  of  0(z),  becomes 

d2<i>/dz2  +  {-%+X/z+(%-m2)/z2}4>  =  0  (6) 

where  6  =  2{£2+K2p2(l-23a2/12)}\  \  =  K2p2a2/6,  p2  h  (l-K2p2a2)/4, 
and  a2  e  r\ 2/(£-Q)2  . 
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Equation  (6)  is  Whittaker's  equation  whose  general  solution  can  be 

expressed  in  terms  of  the  functions  W.  (z)  and  W  .  (-z).  Application 

A,p  -A,|J 

of  the  boundary  and  jump  conditions  yields  the  following  dispersion 
relation: 


W|6WX,M(za)-hWX,M<Za)1+8[llU+6V1  =  ° 


(7) 


“here 

-WX,M(2b)W-X>M(-Za)-  VSWX,M(zb)W-'x,M(-Za)-WX,M(za>K-X>M(-zb)' 


h=xaI^(xa)/IACxa)-2A/fl  and 


2  /,  2> 


g=W^(j(zb)-2it  (zb)/6[£(l-a  /!>  )-Q] 


(zbV 


I£(xb)K£(xd)'I£(xd)Ki(xb) 


The  prime  denotes  differentiation  with  respect  to  the  argument.  The 
subscripts  a,  b,  and  d  denote  evaluation  of  x  and  z  at  r  =  a,  b,  and  d, 
respectively.  A  numerical  evaluation  of  this  dispersion  relation 
exhibits  the  finite  bandwidth  in  k  discussed  earlier,  and  it  is  found 
that  Ak  scales  as  1/rj  and  is  bounded  by  Ak<2 £/r). 

The  dispersion  relation  given  in  Eq.  (7)  is  for  the  beam  frame. 
To  obtain  the  dispersion  relation  in  the  lab  frame,  the  converse  rela¬ 
tion  of  Eqs.  (1)  can  be  used  to  express  Eq.  (7)  in  terms  of  the  lab 

frame  quantities.  In  dimensionless  variables  this  Lorentz  transforma- 
2  -1 

tion  is  Q=y  (Q^-k^  Py  c/u)da)  and  K=y(K£-fi£  py  u^a/c),  where  p=v/c. 
The  effect  of  this  transformation  is  to  shift  the  band  of  unstable 
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modes  to  larger  K  values  and  to  cause  an  asymmetry  of  the  spectrum 
about  the  most  unstable  mode. 

The  most  unstable  mode  in  the  beam  frame  is  the  k=0  mode.  This  is 
true  because  the  damping  term  in  Eq.  (3)  vanishes.  Also,  it  can  be 
argued  that  the  frequency  of  the  k=0  mode  is  either  a  maximum  or  min¬ 
imum  because  in  the  beam  frame  there  is  no  preferred  z-direction  and 
the  dispersion  relation  must,  therefore,  be  an  even  function  of  k. 
Equation  (1)  shows  that  if  the  most  unstable  mode  occurs  at  1^=0,  if 
must  occur  at  k^O.  An  approximate  expression  for  u)  and  k  of  the  most 

unstable  mode  in  the  lab  frame  can  be  found  by  expanding  Eqs.  (1)  about 

2 

k^=0  and  requiring  k^  to  be  real.  Neglecting  terms  of  order  [Im(k^)] 

one  finds  k£=yRe(mb)p/c  and  LO^=yRe(tab)+iIm(wb)/ y.  Because  the  dioco- 

tron  mode  frequencies  are  proportional  to  the  beam  density,  another 

factor  of  1/y  enters;  i.e.,  the  most  unstable  mode  in  the  lab  frame 

2 

occurs  for  k^=Re(u)^)p/c  and  u)^^Re(u)^)+iIm(tu^)/y  where  uj^  denotes  the 

diocotron  mode  frequency  of  Levy's  quadratic  dispersion  relation 

evaluated  using  the  lab  frame  density.  These  equations  suggest  that  as 

2 

far  as  propagation  is  concerned  the  instability  rapidly  (a  1/y  ) 
becomes  less  of  a  problem  as  the  beam  energy  is  increased.  Further¬ 
more,  because  the  scaling  of  the  real  part  of  the  frequency  is  not 
degraded  by  increasing  y,  high  frequency  microwave  generation  using 
intense  high  energy  beams  appears  to  be  possible. 

It  should  be  noted  that  the  relativistic,  finite  k  analysis  pre¬ 
sented  here  expands  the  parameter  space  of  the  diocotron  instability 
problem.  In  addition  to  the  mode  number  &  and  the  geometrical  factors 
b/a  and  d/a,  new  parameters  now  include  K,  y,  r\  and  u^a/c.  The  range 
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of  K  for  unstable  waves,  as  mentioned  earlier,  is  determined  by  r).  The 
parameter  w^a/c,  which  was  introduced  by  requiring  (k,iu))  to  be  a 
four-vector,  is  a  measure  of  how  strongly  the  relativistic  effects 
modify  the  characteristics  of  the  instability.  A  complete  parameter 
study  is  too  involved  to  be  presented  in  a  letter  and  will  be  presented 
elsewhere  along  with  the  results  of  particle-in-cell  simulations  of  the 
diocotron  instability. 

The  authors  wish  to  thank  J.  R.  Cary,  B.  B.  Godfrey,  R.  J.  Faehl, 
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by  the  United  States  Department  of  Energy  and  the  Air  Force  Office  of 
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FIGURE  CAPTION 


Fig.  1.  (a)  Beam  configuration;  (b)  Real  part  of  frequency  vs. 

azimuthal  mode  number  showing  surface  wave  coupling  for 
b/a  =  1.1,  d/a  =  1.5,  k  =  0. 
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